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ABSTRACT 

We have obtained low-dispersion spectra of NGC 3227 with the Space Tele- 
scope Imaging Spectrograph (STIS) to study the intrinsic UV absorption and 
the reddening of the nucleus in this Seyfert 1 galaxy. The UV spectra show a 
wealth of absorption lines at the systemic redshift that span a wide range in 
ionization state (Mg I to N V). The equivalent widths of the lines are consistent 
with our earlier prediction that a "lukewarm absorber" (Te = 18,000 K at the 
ionized face) with a substantial column of gas (N// = 2 x 10^^ cm~^) is present 
and likely responsible for the reddening of the nucleus. The lukewarm absorber 
is also responsible for most of the absorption in the X-rays at energies less than 
1 keV, although a more highly ionized "warm absorber" is needed to account for 
the O VII and O VIII ionization edges. In addition, we require a small column 
(N// = 5 X 10^^ cm~^) of cold gas to match the strengths of the neutral and 
singly-ionized lines in the UV spectra. NGC 3227 is the first Seyfert galaxy in 
which a strong link between the reddening and intrinsic UV absorption has been 
found. By comparing our STIS UV and optical spectra with those of the unred- 
dened Seyfert NGC 4151, we have determined a reddening curve for the nuclear 
continuum source in NGC 3227 over the 1150 - 10,200 A range. The reddening 
curve does not show a 2200 A bump, and is steeper in the UV than reddening 
curves derived for the Galaxy, LMC, and SMC, suggesting a preponderence of 
small dust grains near the nucleus. 
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1. Introduction 

The nearby Seyfert 1.5 galaxy NGC 3227 (z = 0.00386) has been studied extensively 
in most regions of the electromagnetic spectrum. The UV and optical emission lines and 
continuum from the active nucleus are known to be reddened by dust in the Seyfert galaxy, 
although the derived values show a large dispersion (Komossa & Fink 1997a). For example, 
Cohen (1983) determined Eb-v — 0.51 from the narrow Ha/H/? ratio, whereas Winge et al. 
(1995) used the total (narrow + broad) Ha/H/? ratio to obtain Eb-v — 0.28 (the Galactic 
reddening is only Eb-v = 0.02, Schlegel et al. 1998). lUE spectra show that the UV 
continuum is heavily reddened, as evidenced by the steep drop in continuum flux from the 
near to far UV (Courvoisier & Paltani 1992). X-ray observations of NGC 3227 with ROSAT 
and ASCA reveal the presence of an X-ray warm absorber, characterized by an ionization 
parameter (number of ionizing photons with energies > 1 Ryd per hydrogen atom at the 
ionized face) of U = 2.4 and hydrogen (H I plus H II) column density Nh = 3 x 10^^ cm~^ 
(George et al. 1998). lUE observations suggest the presence of intrinsic UV absorption in 
Mg II, although these observations cannot rule out a Galactic origin (Ulrich 1988). 

As pointed out by Komossa & Fink (1997a), the intrinsic neutral column of gas in 
the hne of sight to NGC 3227 (N^/ < 3 x 10^° cm'^ see also George et al. 1998) is not 
sufficient to produce the observed reddening, assuming a Galactic dust/gas ratio (^hi = 5.2 
X 10^^ Eb-v cm~^, ShuU & van Steenberg 1985). Thus, Komossa & Fink suggest that the 
reddening arises in the much higher column of highly ionized gas responsible for the X-ray 
absorption. Similar cases for "dusty warm absorbers" have been made for the QSO IRAS 
13349-1-2438 (Brandt, Fabian, & Pounds 1996) and the Seyfert galaxies NGC 3786 (Komossa 
& Fink 1997b), IRAS 17020+4544 (Leighly et al. 1997), and MCG -6-30-15 (Reynolds 
et al. 1997). 

In a recent paper on the reddening in X-ray absorbed Seyfert 1 Galaxies (Kraemer et 
al. 2000, hereafter K2000), we presented an alternative to the dusty warm absorber model. 
In our model, the dust exists in a "lukewarm absorber" in which the hydrogen is ionized, 
but the columns of O VII and O VIII are negligible. The lukewarm absorber can be placed 
at relatively large distances (> 100 pc) from the nucleus, and therefore account for not only 
the reddening of the continuum source and broad- line region (BLR), but the reddening of 
the narrow-hne region (NLR) as well. 

In K2000, we showed a specific model for a lukewarm absorber in NGC 3227. To 
constrain the model, we required 1) the distance from the continuum source to be greater 
than ~100 pc, to account for the reddening of the majority of the NLR emission, and 2) the 
column density to be N^^ = 2 x 10^^ cm~^, to account for the average reddening of Eb-v 
— OA from the hterature (assuming the Galactic dust to gas ratio). The parameters of the 
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lukewarm absorber were then determined from a fit to the soft X-ray absorption observed 
by ASCA and ROSAT (see K2000 for details). The resulting model was characterized by 
the parameters U = 0.13, Uh = 20 cm~^, and distance = 120 pc. Interestingly, most of 
the soft X-ray absorption can then be attributed to He II Lyman continuum absorption. 
Our lukewarm absorber model for NGC 3227 predicted large ionic columns of N V, C IV, 
Si IV, and Mg II, which we suggested could be tested with observations in the UV at higher 
sensitivity than possible with lUE. 

As part of an ongoing investigation into the intrinsic UV absorption in Seyfert galaxies, 
wc obtained spectra of NGC 3227 with the Space Telescope Imaging Spectrograph (STIS) 
on the Hubble Space Telescope {HST) under a STIS guaranteed time observations (GTO) 
program. Given the HST time available for this object (2 orbits), the UV spectrum is too 
faint to obtain reasonable signal-to-noise levels with the echelle gratings, so we opted for 
low-dispersion UV spectra through a long sht (a description of STIS and its instrumental 
configurations can be found in Lcitherer et al. 2000). We also obtained optical spectra 
through the same slit to investigate the reddening of the continuum over a broad wavelength 
region. 



2. Observations and Measurements 

We observed the nucleus of NGC 3227 on 2000 February 8 using a 52" x 0".2 slit and 
the STIS low-dispersion gratings. Wc show the details of the observations in Table 1. We 
reduced the spectra using the IDL software developed at NASA's Goddard Space Flight 
Center for the STIS Instrument Definition Team. The nucleus dominates the light in the 
long-slit images. We extracted the nuclear spectra using a 0".3 height (perpendicular to the 
dispersion). We then combined the reduced spectra from each grating in their regions of 
overlap to produce a final spectrum over the 1150 - 10,200 A range at a spectral resolving 
power of A/AA 1000. 

Figure 1 shows the STIS UV/optical spectrum of the nucleus in NGC 3227 and that of 
the Seyfert 1.5 galaxy NGC 4151 (Nelson et al. 2000), scaled by a factor of 0.15 in fiux. The 
severe extinction in the UV spectrum of NGC 3227, compared to NGC 4151, is immediately 
obvious. The smooth curves in Figure 1 are spline fits to continuum regions that are free of 
any obvious contamination by absorption or emission lines including the "little blue bump" 
between 2200 and 4200 A in the NGC 4151 spectrum, which consists of a blend of broad 
Fe II and Balmer recombination emission (Wills, Netzcr. & Wills 1985). This feature may 
also be present in the NGC 3227 spectrum, but since it cannot be distinguished as a separate 
feature from the continuum emission, the continuum fit in this region may be considered an 
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upper limit. 

In Figure 2, we show an expanded view of the UV spectrum of NGC 3227. In addition 
to the broad and narrow emission hnes and the interstellar absorption lines from the Galaxy, 
a number of intrinsic absorption lines are present and identified. These absorption lines 
appear at the systemic redshift of NGC 3227 (z = 0.00386 from H I 21 cm observations 

of the host galaxy, dc Vaucouleurs et al. 1991), and span a wide range in the ionization 
potential needed to create these ions, from zero (O I, C I, Mg I) to 77.5 eV (N V). All of 
the absorption lines are unresolved at the velocity resolution of the spectra (~ 300 km s~^ 
FWHM). 

To measure the absorption lines, we fit cubic splines to adjacent wavelength regions, 
which typically contain both continuum and emission-line contributions. Since the absorp- 
tion lines are unresolved, the column densities cannot be determined directly from integration 
of optical depths (e.g., Crenshaw et al. 1999), and we must therefore rely on the equiva- 
lent widths and radial velocity centroids (relative to the systemic redshift). A few of the 
absorption lines that we have identified are blended with stronger lines, and we are not able 
to determine rehable values for these hnes. For the blended C IV doublet, we were able 
to model the observed blend with two lines of approximately equal equivalent width. We 
show our absorption-line measurements in Table 2 (along with column densities from models 
discussed in section 4). 



3. The Reddening Curve for NGC 3227 

To determine a reddening curve from the nuclear continuum emission in NGC 3227, we 
assumed that 1) the intrinsic UV/optical continuum of NGC 4151 is essentially unreddened, 
and 2) the relative flux of NGC 4151 as a function of wavelength represents the intrinsic 
continuum spectrum of NGC 3227. The flrst assumption is based on the finding by Kriss et 
al. (1995) that the far-UV continuum of NGC 4151 is consistent with a power-law plus a 
small amount of reddening that is due solely to our Galaxy {Eb-v = 0.04, Burstein & Heiles 
1982). The second assumption is based on evidence that the intrinsic continua of Seyfert 
1 galaxies have a relatively standard shape, which is modified by the effects of dust and 
contributions from the host galaxy to produce the observed spectral energy distributions 
(Ward et al. 1987). We will explore the uncertainties associated with our assumptions later 
in this section. 

Given our assumptions, we can determine the value of E(B— V) = 2.5 [log(XB) — log(Xy)], 
where X is the ratio of the continuum fit to NGC 4151 to that of NGC 3227 as a function of 
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wavelength (Figure 1), and is evaluated at the effective wavelengths of the B (4400 A) and V 
(5500 A) filters. Thus, we find that E(B-V) = 0.18. We can also determine the reddening 
curve at any wavelength, relative to V and as per convention normalized to E(B— V) (Savage 
& Mathis 1979): 

E{X - V) ^ A, -Ay ^ logjXx) - logjXy) 

E{B -V) Ab~ Ay logiXs) - log{Xy) ' ^ ^ 

where A^, Ay, and A^ are the extinctions in magnitudes at B, V, and an arbitrary A. 

Since we do not know the intrinsic continuum flux of NGC 3227 at any wavelength, 
we cannot determine the total extinction directly. Thus, we need to calibrate the above 
reddening curve to determine the actual extinction law: 

(2) 



E(B - V) E(B - V) 

by determining the offset Ry = Ay/E(B— V). The value of Ry for the "standard" 
Galactic reddening curve is 3.1 (Savage & Mathis 1979), but this value can vary among 
different lines of sight in the Galaxy (Cardelli, Clayton, & Mathis 1989). However, extinction 
laws appear to have essentially the same values at wavelengths A > 7000 A (Cardelli et al. 
1989). If we assume this is the case for NGC 3227, then Ry is the constant we need to add 
to E(A— V)/E(B— V) to match the standard Galactic extinction law at long wavelengths. 

Figure 3 shows extinction laws for the Galaxy, LMC, and SMC, along with the extinction 
law determined for NGC 3227 in the fashion described above. The shape of the extinction 
curve for NGC 3227 is very similar to the others over the 4000 - 10000 A range, and we find 
that an offset of Ry = 3.2 provides a good match to these curves; this value is very close to 
the standard Galactic value of Ry =3.1. 

The extinction curves in Figure 3 begin to differ strongly at A < 4000 A. The standard 
Galactic curve (Savage & Mathis 1979; Seaton 1979) shows a strong 2200 A bump, and 
the weakest (although substantial) rise to the UV. The LMC curve (Koornneef & Code 
1981; Fitzpatrick 1985) shows a stronger rise to the UV and a less pronounced bump. The 
SMC curve shows the strongest rise to the UV (until now) and no detectable 2200 A bump 
(Hutchings et al. 1982; Witt & Gordon 1999). 

The extinction curve for NGC 3227 is extreme, in that it shows a steeper rise to the UV 
than all of the other curves. Furthermore, the extinction curve resembles that of the SMC, 
in that there is no 2200 A bump. Interestingly, Pitman, Clayton, & Gordon (2000) find a 
similar result for higher luminosity AGN; specifically, there is no convincing evidence for a 
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2200 A feature in the spectrum of any QSO. Tests show that lower continuum fits in the 
2200 - 4200 A region (assuming a substantial little blue bump, see section 2 and Figure 1) 
have little effect on the overall sharp rise to the far-UV. Given the reddening of E(B— V) = 
0.18, we can determine the extinction of NGC 3227 at any wavelength from Figure 1: e.g., 
Aa(1200 A) = 5.3, Aa(2000 A) = 2.9, Aa(3000 A) = 1.2, and Ay(5500 A) = 0.58. With this 
extinction curve, we find that the intrinsic continuum flux of NGC 3227 at any wavelength is 
then 0.17 times that of NGC 4151 at the respective times of observation. Note, however, that 
the STIS spectrum of NGC 3227 was obtained at a relatively low continuum state (section 
4.2), whereas that of NGC 4151 was obtained at a moderately high state (Crenshaw et al. 
2000). 

To test the uniqueness of our reddening curve for NGC 3227, we have explored possi- 
ble systematic uncertainties associated with our assumptions. The assumption that NGC 
4151 represesents the intrinsic continuum spectrum of NGC 3227 is based on the observed 
similarity of UV/optical continua in Seyfert 1 galaxies (e.g.. Wills, Netzer, & Wills 1985). 
Unfortunately, HST observations of other unreddened Seyfert 1 galaxies covering the full 
UV/optical range at high spatial resolution (to avoid galaxy contamination) are not yet 
available. Therefore, we use the extreme variations of NGC 4151 itself, which was monitored 
extensively with lUE. Clavel et al. (1990) find that the ratio of UV to optical continuum 
Fa(1450) /Fa(5000) varies by a factor of 3.1, with high ratios corresponding to high UV fluxes. 
However, Kaspi et al. (1996) find that a portion of this effect is due to a substantial galaxy 
contribution in the apertures used for optical observations; subtracting the stellar fiux at 
5000 A leads to a ratio that varies by a factor of 2.1. Thus, if we assume an extremely high 
state for NGC 4151 and an extremely low state for NGC 3227, this factor would bring the 
reddening curve of NGC 3227 close to (but still above that) of the SMC, since the reddening 
curves of these two objects in Figure 3 differ by 1.0 mag (on average) in the 1200 - 2000 A 
range (for E^-v = 0.18). Other systematic uncertainties will tend to steepen the reddening 
curve slightly. Assuming that the entire 2500 - 4000 bump in the observed spectrum (Figure 
1) is due to the little blue bump leads to an artificial bump in the extinction curve that spans 
this region (with an amplitude of 0.3 mag), which does not resemble the narrower 2200 A 
bump in the standard Galactic curve. Scattering of the radiation back into our line of sight 
could potentially "bluen" the continuum, and correction for this effect would steepen the 
reddening curve in the UV. However this effect must be small, since our extraction slit covers 
only a 15 pc X 22 pc region, and we have estimated the dusty gas cloud is ~100 pc from the 
nucleus. Assuming a spherical shell with a covering factor of one and uniform scattering, 
only ~ 1% of the UV light scattered in our direction enters the aperture. Furthermore, dust 
grain models show that most of the UV radiation is absorbed and re-emitted in the infrared, 
rather than scattered (Draine & Lee 1984). Thus, our estimate of the systematic uncertain- 
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ties show that we cannot rule out a reddening curve for NGC 3227 that is essentially the 
same as that of the SMC, although it is hkely to be steeper. 

The STIS spectra also provide an opportunity to determine the reddening of the emission 
lines from the inner NLR (within 0".15 or 11 pc of the nucleus for Hq = 75 km s~^ Mpc~^) 
for comparison with the continuum reddening. At the temperatures and densities of the 

NLR, the He II A1640 and A4686 lines are due to recombination, and the intrinsic He II 
A1640/A4686 ratio is ~7.2 (Seaton 1978; Kraemer et al. 1994). We were able to deblend 
the narrow components of these lines from their broad counterparts, using the [O III]A5007 
lines as a template, to determine an observed ratio of 0.19 ±0.05. This yields Aa(1640) — 
A;^(4686) = 3.9 ±0.3 mag for the narrow lines. By comparison, the continuum extinction 
curve yields A;^(1640) — A;v(4686) = 3.4. These values are in reasonably good agreement, and 
provide further evidence that the reddening of the continuum and emission lines is primarily 
due to material outside of the inner NLR (i.e., our lukewarm absorber). The small additional 
reddening experienced by the emission lines may be due to dust in the inner NLR, which 
has been detected in other Seyfert galaxies (e.g., NGC 1068, Kraemer & Crenshaw 2000). 

Our value for E(B— V) in NGC 3227 (= 0.18) is somewhat lower than previous estimates 
(section 1). However, previous estimates have assumed the standard Galactic extinction 
curve, which is clearly incorrect, and the E(B— V) values are therefore sensitive to the wave- 
lengths of the reddening indicators Tised. Also, the use of narrow Ha/H/? is complicated by 
the deconvolution of the narrow and broad components, and in particular, the deconvolution 
of narrow Ha and [N II] AA 6548, 6884, which may account for at least part of the large 
dispersion in E(B— V) values. (Unfortunately, we cannot obtain a narrow Hq:/H/3 ratio from 
our own observations, because the spectral resolution of our own data does not allow an 
accurate deconvolution of the Ho; and [N II] lines.) Finally, we note that previous reddening 
measurements were obtained for much larger projected regions (> 1"), and the reddening 
may be somewhat aperture dependent. 

4. The Intrinsic UV Absorption Lines 

4.1. Properties of the Absorbers 

The high-ionization absorption lines in Figure 2 show significant saturation, suggesting a 
large column of gas. In particular, the C IV AA1548.2, 1550.8 lines arc blended and approach 
zero intensity in their troughs, despite the low resolution which decreases the observed depths 
of the cores of unresolved absorption hues. The N V AA1238.8, 1242.8 and S IV AA1393.8, 
1402.8 absorption lines in Figure 2 also suggest a substantial column of ionized gas, in 
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that the longer wavelength member of each doublet is much stronger than the 1:2 ratio for 
unsaturated lines. Qualitatively, these lines confirm our prediction in K2000 that a lukewarm 
absorber is present in NGC 3227. 

Low-ionization absorption lines are relatively rare in Seyfert 1 galaxies; in an HST 
survey, we found that only one of 10 objects with intrinsic C IV and N V absorption showed 
Mg II in absorption (Crenshaw ct al. 1999). This object is NGC 4151, which shows not 
only low-ionization species, but absorption from fine-structure and metastable levels as well 
in most kinematic components (Kraemer et al. 2001, and references therein), indicating a 
range of moderate to high densities (n,. from 10^ cm^"^ to >10^ cm~^). For NGC 3227, the 
fine structure Unes of O I, C II, Si II, and Fe II are not detected, indicating low densities (ne 
< 100 cm~^). Note that our original model in K2000 predicted negligible ionic columns of 
C I and Fe II. Our detection of these fines suggests the presence of a small column of cold 
gas, as discussed in subsequent sections. 

The intrinsic absorption fines listed in Table 2 are near the systemic redshift of the host 
galaxy: the average radial velocity of all of the lines is -|-20 ± 125 km s~^. Although most 
intrinsic absorbers in Seyfert 1 galaxies are in outflow, systems near the systemic redshift 
are not uncommon (Crenshaw et al. 1999). There is some evidence that the high-ionization 
lines are at a sligfitly different radial velocity tfian tfie low-ionization lines: the average radial 
velocity for N V, C IV, and Si IV is —70 ± 55 km s""*^, whereas the average radial velocity 
for the remaining, low-ionization lines is -1-125 ± 50 km s^^. Future observations at fiigfier 
dispersions would be helpful in identifying multiple kinematic components and confirming 
this velocity difference. 

4.2. Photoionization Models of the Absorbers 

We have generated photoionization models to match the observed equivalent widths of 
the absorption lines. The photoionization code has been described in previous publications 
(cf. Kraemer et al. 1994). The photoionization models are parameterized in terms of the 
number of ionizing photons per hydrogen atom at the ionized face of the absorber ([/), the 
density (n//), the distance from the ionizing source (D), and the sum of the H I and H II 
column densities (N/f). We assume that gas is ionized by the continuum radiation emitted 
by the central source in the active nucleus of NGC 3227. As discussed in K2000, we modeled 
the EUV to X-ray spectral energy distribution (SED) as a series of power-laws of the form F,^ 
oc z^~", with a = 1 below 13.6 eV, a = 2 over the range 13.6 eV < hu < 500 eV, and a = 0.6 
above 500 eV. Note that the spectral index above 500 eV matches the value derived from 
the best fit to the combined ROSAT/PSPC and ASCA datasets (George et al. 1998). Based 
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on this SED, and the assumption that NGC 3227 and NGC 4151 have similar optical-UV 
SEDs, we derived a luminosity in ionizing photons, from 13.6 eV to 10,000 eV, of ~ 1.5 x 
lO'"^''' photons s~^. We assumed roughly solar element abundances (cf. Grevesse & Anders 
1989), which are, by number relative to H, as follows: He = 0.1, C = 3.4 x lO^'^, N = 1.2 x 
10-4, O = 6.8 X 10-^ Ne = 1.1 x 10~^ Mg = 3.3 x 10"^ Si = 3.1 x 10-^ S = 1.5 x 10-^ 
and Fe = 4.0 x 10"^. Further, we assumed that both silicate and carbon dust grains are 
present in the gas, with a power-law distribution in sizes (see Mathis, Rumpl, & Nordsieck 
1977). For our original lukewarm model, we assumed depletions of elements from gas phase 
typical of the Galactic interstellar medium (cf. Snow & Witt 1996): C, 65%; O, 50%; Si 
95%; Mg 82%; and Fe 95% (K2000). 

In K2000, we determined that the observed reddening and the bulk of the X-ray ab- 
sorption could be explained by a single component of ionized gas (a lukewarm absorber), 
which covers the continuum source, broad-line region, and much of the narrow-line region. 
As discussed in Section 1 and K2000, we had previously fixed the total column density of the 
lukewarm absorber at N// = 2 x 10^^ cm"^, based on the reddening derived from the Balmer 
decrement and the Galactic dust-to-gas ratio. From our empirical reddening curve for NGC 
3227, we have determined E(B— V) = 0.18, roughly half the value we assumed in K2000. 
and consequently one might argue that Nj^ should be reduced by the same factor. However, 
since the reddening curve is quite different in NGG 3227 than in the Galaxy, it is unlikely 
that the Galactic relationship of the reddening to the total hydrogen column strictly applies. 
Furthermore, our success in fitting the soft X-ray absorption gives us confidence that the 
He H column density, and hence the U and N^^ are roughly correct. Therefore, we have 
opted to keep the same parameters as before for the lukewarm component: U = 0.13, uh — 
20 cm~^, and N// = 2 x 10^^ cm^^. Note that the distance of the lukewarm absorber from 
the central source depends on the luminosity, and we derived D = 120 pc from our original 
finding that the reddening-corrected luminosity of NGC 3227 derived from the 1992 spectra 
in Winge et al. (1995) is equal to that of NGC 4151 from Nelson et al. (2000). However, 
the smaller reddening derived in this paper places the luminosity of NGG 3227 at 0.59 times 
that of NGG 4151, which yields D = 92 pc. Furthermore, we have assumed that the 1992 
luminosity represents the time-averaged value, although the current STIS spectrum indicates 
that it is in a lower continuum state than in previous observations (i.e., the continuum flux 
at 5500 A is 30% that of the average value from the 1992 Winge et al. spectra). 

As noted in section 4.1, the STIS spectra reveal the presence of saturated N V and 
G IV lines, in qualitative agreement with the predictions of the lukewarm model. However, 
the equivalent widths of the Si IV AA1393, 1403 lines arc about twice those expected from 
our original model and a curve-of-growth that fits the C IV and N V lines (see the next 
subsection), which is evidence for a greater column of Si IV than we had predicted. We 
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suggest that this is the result of a smaller fractional depletion of silicon onto dust grains. 
This is not surprising, given the unique extinction law we have determined for NGC 3227. 
However, since we do not have an independent way to determine the elemental depletions 
for the constituents of the silicate grains, we have recomputed the lukewarm model using 
the Galactic depletions, scaled by 0.8 to account for the larger fraction of silicon in gas- 
phase, as follows: O, 41%; Si, 77%; Mg, 67%; and Fe, 78%. While the 2200 A bump is 
generally assumed to be due to carbon-based grains, the origin of this feature is not well 
understood (see Duley & Seahra 1999, and references therein) and its relative strength may 
not be correlated with the amount of carbon depletion (Cardelli et al. 1996). Hence, we 
have assumed the same depletion of carbon as in K2000. 

As we will demonstrate in section 4.3, the lukewarm model predicts relatively small ionic 
columns for most low- ionization species (and essentially no Fe II), although such lines are 
strong in our STIS spectra (Table 2). We attempted to increase the low ionization column 
densities by increasing Nj^, but models that fit the observed EW's of the Si II, Mg II, 
and Fe II lines underpredict those of the O I and C II lines. Furthermore, there is some 
evidence that the low ionization lines have a different velocity radial centroid than that of the 
high ionization lines (section 4.1). Hence, we have included an additional component of low 
ionization gas (a "cold component" ). Based on the absence of fine structure lines of Si II and 
C II, the density of the cold component is < 100 cm~^ and, consequently, cannot be closer to 
the source of ionizing radiation than the lukewarm gas (for example, if the two components 
were at the same radial distance, the cold component would still have a fairly high ionization 
parameter of U = 0.03). Therefore, we have placed the cold component further from the 
nucleus and, as a result, it is screened by the lukewarm gas. In modeling this component, 
we have assumed the same abundances and depletions as the lukewarm component, and the 
following parameters: U = 0.00015, tlh = 100 cm^^, = 5.2 x 10^^ cm^^, and D = 160 pc. 
The column density was set to provide the best fit to the observed low ionization lines, as 
described in the next section. The resulting ionic column densities for both lukewarm and 
cold components are given in Table 2. 



4.3. Comparison of Observations and Models 

To demonstrate that the equivalent widths of the absorption lines are consistent with 
the predictions of our absorber models, we have generated Maxwcllian curves of growth for 
different velocity widths ("b" values in km s^^, see Spitzer et al. 1978), as shown in Figure 
4. In this figure we overplot the absorption-line data from Table 2 (in a few cases, lines 
from the same ion have nearly identical values and only one line is plotted for clarity). For 
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each absorption line that is overplotted, the vertical position is determined by the observed 
equivalent width and the horizontal position is determined by the ionic column predicted by 
the model. For a group of lines to be consistent with a model, the lines should lie relatively 
near a single curve of growth representing a particular b value. 

The top panel in Figure 4 shows our results for the lukewarm absorber model described 
in the previous section. This plot shows that the high-ionization lines (C IV, N V, and Si IV) 
arc consistent with the ionic column densities from the lukewarm model and a b value of 
100 ± 50 km s^^. Mg II also matches the curve of growth at b = 100 km s^^. However, the 
column densities of the low-ionization lines of C II, Si II and Mg I are underprcdictcd, and 
the lukewarm model predicts negligible columns of O I, C I, and Fe II. Thus, an additional 
component of relatively cold gas is required to match these lines. 

The bottom panel in Figure 4 shows our result for the cold absorber model described 
in the previous section. Nearly all of the low-ionization lines are consistent with the ionic 
columns predicted by the model and b = 200 ± 50 km s~^. However, the columns of Mg I and 
Mg II are overpredicted, which is exacerbated by the fact that some of the Mg II equivalent 
width must be attributed to the lukewarm absorber. This suggests that the depletion of 
Mg into dust may actually be somewhat larger than the value chosen for our models (i.e., 
85% rather than the 67% we assumed), and that, in general, the relative depletion of the 
constituents of the dust grains differs from that in the Galactic ISM. The only low-ionization 
line that is severely underprcdictcd by this model is C I. Our underprediction of the C I 
column is likely the result of the unusually steep extinction curve in the far UV, which is 
not accounted for in our model of the lukewarm absorber. The large extinction inferred 
at ~1100 A depletes photons with energies near the ionization potential of C I (11.26 eV), 
which would likely result in a substantial C I column within the slab (note only ~1% of the 
C needs to be C I.) 

Overall, the curves of growth in Figure 4 show that the equivalent widths of the absorp- 
tion lines are consistent with our models of the lukewarm and cold absorbers in the line of 
sight to the nucleus. The different b values and possibly different velocity centroids for the 
low and high-ionization lines are further evidence for two components with different kine- 
matics. Spectral observations of the absorption at higher dispersion, although challenging 
due to the low UV flux, would be helpful in providing further tests of and/or refinements to 
these models. 



5. Conclusions 



We have shown that there are three very different intrinsic components of gas in the 
hne of sight to the nucleus of NGC 3227: an X-ray warm absorber (U = 2.4, Nh = 3 x 10^^ 
cm^^), a lukewarm absorber (U = 0.13, Nh =2 x 10^^ cm^^), and a cold absorber (U = 1.5 
X 10"'*, = 5 X 10^^ cm~^. The X-ray warm absorber, characterized by O VII and O VIII 
absorption edges, shows a variable column of highly ionized gas, as discussed in detail by 
George et al. (1998). The lukewarm absorber predicted by K2000 and revealed by the STIS 
spectra shows high-ionization lines (N V, C IV, Si IV) typical of intrinsic UV absorbers 
in most Seyfert 1 galaxies (Crenshaw et al. 1999). The low ionization lines in the STIS 
spectra are primarily due to the cold component, although our model indicates significant 
contributions to the observed C II, Si II, and Mg II from the lukewarm absorber. The H I 
column from the cold component is only 3 x 10^^ cm~^ (Table 2), which is much smaller 
than the intrinsic neutral column of ~3 x 10^° cm~^ derived from the soft X-ray absorption 
(Komossa & Fink 1997a; George et al. 1998). However, as we pointed out in K2000, most of 
the soft X-ray absorption can be attributed to He II Lyman continuum absorption from the 
lukewarm absorber, and thus previous estimates of the neutral column from X-ray studies 
of NGC 3227 are hkely a factor of ten too large. 

We find that our lukewarm absorber has sufficient column to provide the dust required 
to redden the continuum and emission lines in NGC 3227, assuming a Galactic dust-to- 
gas ratio. Furthermore, the lukewarm absorber is at a distance sufficient to account for the 
reddening of the NLR. Thus, we prefer the lukewarm absorber as the source of the reddening. 
Although a dusty warm absorber cannot be completely ruled out, the variability of the warm 
absorber in NGC 3227 over a time scale of two years suggests a location that is within light 
days of the central source (George et al. 1998), which is unsuitable for reddening the NLR. 
Further monitoring would be helpful in verifying the variability and establishing the location 
of the X-ray absorber. 

We note that additional support for reddening by dust at large distances comes from 
spectropolarimetry of NGC 3227. Schmidt & Miller (1985) find that the continuum polariza- 
tion increases from 1% to 3% going from 6000 A to 3500 A. Furthermore, narrow emission 
lines are polarized by the same or slightly smaller amounts as nearby continuum regions, 
suggesting scattering by ahgned dust grains at the location or outside of the NLR. This 
confirms the suggestion by Thompson et al. (1980) that the nucleus, BLR, and NLR of 
NGC 3227 are covered by a polarizing dust cloud. 

The extinction curve that we have derived for NGC 3227 has important imphcations for 
reddening corrections to Seyfert galaxies. For example, photoionization models of the NLR 
rely heavily on UV/optical fines ratios (e.g., Lya/H/?, C IV/H/3) to determine the physical 
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conditions in the emitting gas (cf., Kraemer & Crenshaw 2000). However, the determination 
of intrinsic hne ratios is very sensitive to the reddening curve, which is almost always assumed 
to be Galactic. For E^_y = 0.18 and the extinction curve for NGC 3227 in Figure 3, the 
application of the standard Galactic curve would result in an underestimate of the extinction 
at 1200A of 3.6 mag (i.e., the corrected flux would be too small be a factor of 27.5). Note 
that this effect is mitigated if the He II A1640 and He II A4686 lines are used to determine 
the reddening; the long UV/optical baseline provides a crude first-order correction, even if 
the wrong reddening curve is assumed. However, the use of the Galactic curve would still 
result in substantial errors, and in particular, an overestimate of the extinction in the 2200 
A bump region. 

The extinction curve also has important implications for the properties of the dust near 
the nucleus, since the shapes of the reddening curves depend on the composition and size 
distribution of the dust grains. Theoretical work in this area (Mathis, Rumpl, & Nordsieck 
1977; Draine & Lee 1984) attributes the ratio of UV to optical absorption to the sizes of 
the grains (large UV absorption corresponds to smaller average grain size). Thus, the large 
UV extinction suggests a smaller average size of dust grains (compared to the Galaxy) for 
the SMC (Pei 1992; Weingartner & Draine 2000) and, by analogy, NGC 3227. One possible 
explanation for the smaller grain sizes near the nucleus of this Seyfert galaxy is that they 
are due to shocks, which tend to destroy large grains via shattering (Jones et al. 1996). 

It is interesting that NGC 4151 harbors a UV absorber with a large column but shows 
no evidence for dust in the gas. STIS cchcllc observations of NGC 4151 show a kinematic 
component of absorption at —500 km s^^ (relative to systemic) that is characterized by U 
= 0.015 - 0.060 (variable) and N// = 2.8 x 10^^ cm-^ (Kraemer et al. 2001), but, as we 
have discussed, the UV/optical continuum of NGC 4151 shows no significant reddening. 
However, the high-column absorber in NGC 4151 is at a distance of only ~ 0.03 pc from 
the nucleus (Kraemer et al. 2001), which is the approximate dust sublimation radius for the 
UV luminosity of this Seyfert (Barvainis 1987). Furthermore, given the large difference in 
distance from the nuclear continuum source, it is possible that the high-column UV absorbers 
in these two Seyferts have completely different origins. 

iif5T observations of the (faint) UV emission from other reddened AGN would be helpful 
in exploring the relationship between absorption and reddening in these objects. Given that 
a substantial fraction of Seyfert 1 galaxies show evidence for reddening of their continua by 
dust (Ward et al. 1997), it would be interesting to determine their reddening curves with the 
procedures we have presented in this paper. New observations would be particularly useful 
for determining the range of intrinsic spectrum shapes in unreddened AGN and exploring 
the variety of extinction laws in reddened AGN. 
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Fig. 1.— STIS UV and optical spectra of NGC 3227 and NGC 4151 (scaled in flux by a 
factor of 0.15). The smooth lines are spline fits to continuum regions in each spectrum. 

Fig. 2. — Expanded view of the UV spectra of NGC 3227. Intrinsic absorption hnes detected 
near the systemic redshift (z = 0.00386) are labeled. 

Fig. 3. — Reddening curve determined for NGC 3227 as a function of wavelength. Reddening 
curves for the Galaxy (the standard curve of Savage & Mathis 1979), LMC (Koornneef & 
Code 1981) and SMC (Hutchings 1982) are given for comparison. 

Fig. 4. — Empirical curves of growth, based on a Maxwellian velocity distribution, for 
different velocity widths ("b" values given in km s~^, Spitzer et al. 1978). The vertical 
position of each point is determined from the observed equivalent width of an absorption 
line. The horizontal direction is determined by the column density for the corresponding 
ion from the model. The upper and lower panels give values based on our models of the 
lukewarm and cold absorption components, respectively. 
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Table 1. STIS Observations of NGC 3227 



Grating 


Exp. Time 


Coverage 




(sec) 


(A) 


G140L 


2132 


1150 - 1715 


G230L 


1612 


1573 - 3157 


G430L 


120 


2900 - 5710 


G750L 


120 


5267 - 10,261 
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Table 2. Intrinsic Absorption - Measurements and Models 



Line 


V " 


EW^ N (Lukewarm)^''^ 




(Inmrnprif^ on 






(A) 


(cm ^) 




Ob^prv^^f lonc; 


Si III A1206.5 






6.8 X lO^'' 


9.4 X 10^2 


blended with Lya 


Lja A1215.7 






1.8 X 10^^ 


3.0 X 10^9 


saturated 


N V A1238.8 


_44 


1.46 ±0.24 


4.8 X 10^^ 






N V A1242.8 


-71 


1.23 ±0.15 


4.8 X 10^^ 






Si II A1260.4 


+ 103 


1.78 ±0.23 


7.7 X 10^3 


3.5 X 10^^ 




C I A1277.2 


+145 


1.07 ±0.39 




3.4 X 10^3 




I A1302.2 


-236 


2.94 ±0.52 


1.1 X 10^^ 


1.3 X 10^^ 


includes Si II A1304.4 


Si II A1304.4 






7.7 X 10^^ 


3.5 X 10^^ 


blended with I 


C I A1328.8 


+174 


0.99 ±0.23 




3.4 X 10^^ 


includes Gal. C II 


C II A1334.5 


+100 


3.22 ±0.24 


1.5 X 10^^ 


5.8 X 10^5 




Si IV A1393.8 


-38 


1.92 ±0.45 


1.1 X 10^^ 






Si IV A1402.8 


-46 


1.43 ±0.20 


1.1 X 10^5 






Si II A1526.7 


-77 


1.12 ±0.19 


7.7 X 10^3 


3.5 X 10^^ 




C IV A1548.2 


-170 


2.38 ±0.24 


7.8 X 10^6 




EW from fit to doublet 


C IV A1550.8 


-170 


2.38 ±0.24 


7.8 X 10^^ 




EW from fit to doublet 


C I A1560.3 


+72 


1.11 ±0.35 




3.4 X 10^3 




C I A1656.9 


-162 


9 1 Q _i_n 




3.4 X 10^3 




Al II A1670.8 


-9 


0.85 ±0.23 








Fe II A2344.2 


+198 


1.89 ±0.27 




4.0 X 10^^ 




Fe II A2374.4 








4.0 X 10^^ 


blended with Gal. Fe II 


Fe II A2382.8 


+151 


2.40 ±0.41 




4.0 X 10 




Fe II A2586.6 








4.0 X 10^^ 


blended with Gal. Fe II 


Fe II A2600.2 


+122 


2.24 ±0.38 




4.0 X 10^^ 




Mg II A2796.3 


+11 


2.66 ±0.32 


1.4 X 10^^ 


3.2 X 10^^ 




Mg II A2803.5 


+115 


1.52 ±0.20 


1.4 X 10^^ 


3.2 X 10^^ 




Mg I A2853.0 


+61 


0.62 ±0.11 


2.5 X 10^2 


1.1 X 10^3 




^relative to the systemic redshift of 0.00386 






''blank entry: blended with another line. 








■^blank entry: not predicted by model. — 


- entry: negligible 


ionic column 


= 0.13, N(H I + H II) 


= 2 X 10^1 cm 


-2 






= 1.5 X 10" 


N(H I + H II) = 5.2 X 


10^9 cm-2. 
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Fig. 1. 



- 21 - 



0.50 



-I— I— I— I— I— I— I- 



o< 



0.40 



0.30 ~- 



m 

CP 

CD 



O 



X 




1200 



1300 



1400 



1500 



Wavelength (A) 



1600 



1700 



o< 



u 



in 

CD 



O 



X 
D 




1800 2000 2200 2400 g600 2800 3000 

Wavelength (A) 
Fie. 2. 



A,/E(B-V) 




Fig. 3. 



-23- 




Model log NfA (cm"') 



Fie. 4. 



